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1. Introduction 

Already the classical works of nineteenth-century geometers established a major 
connection between differential geometry and the theory of partial differential 
equations. Powerful solution-generating techniques such as the Backlund and 
Darboux transformations [33] have origins in the prototypical relationship between 
pseudospherical surfaces and solutions of the sine-Gordon equation. 

Methods available for solving nonlinear partial differential equations were 
substantially extended in the seventies to include the inverse scattering transform and 
its numerous developments; see, e.g., [8, 14, 26, 39]. An important open problem is to 
describe the class of partial differential equations solvable by these powerful methods. 
Indirect detectors such as the symmetry analysis have been involved in obtaining 
extensive complete classifications of integrable evolution equations and systems; see [28] 
and references therein. The known theoretical answer given in terms of the existence of 
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the associated one-parametric zero curvature representation 
Ay-B x + [A, B} = 

has been considered as a classification tool in conjunction with the gauge cohomology 
by one of us [25]. These methods are not limited to evolution equations, although 
the necessary computations are rather complex, resource consuming, and unthinkable 
without substantial use of computer algebra. However, certain partial differential 
equations of geometric origin are particularly well suited for this classification method, 
namely, the Gauss-Mainardi-Codazzi equations of immersed surfaces. These equations 
always possess an associated linear zero curvature representation, albeit without the 
spectral parameter. 

Since their introduction by Weingarten [36], immersed surfaces in R 3 that satisfy a 
functional relation between the principal curvatures have been of continuing interest in 
differential geometry, see, e.g., [18, 20, 22]. It is therefore not surprising that attempts 
have been made to identify classes of Weingarten surfaces such that the corresponding 
Gauss equation is integrable in the sense of soliton theory. The work of Wu [38] and 
Finkel [16] indicated that all integrable cases are classical, characterized by a linear 
relation between the Gauss and the mean curvatures (linear Weingarten surfaces [12, 
§812]; see also [17] and references therein). In other words, the integrable Weingarten 
surfaces were conjectured to be either minimal or parallel to surfaces of constant 
Gaussian curvature. This conjecture was, however, disproved by the present authors 
in [1], henceforth referred to as Part I. In Part I we found another integrable class, 
consisting of surfaces with a constant difference between the principal radii of curvature, 
which we called surfaces of constant astigmatism. Surprisingly enough, this extra class 
turned out to be classical as well, apparently first mentioned by Beltrami [3, Ch. 9, §20], 
covered by Bianchi [4] and Darboux [12], see also [31], yet forgotten today. 

In this paper we continue the work begun in Part I and complete the classification of 
integrable classes in the simplest possible case. The integrability criterion we adopt is the 
existence of an s 1(2)- valued zero curvature representation depending on a nonremovable 
parameter. We apply the same method of formal spectral parameter, introduced in [25] 
and briefly reproduced in Part I. The underlying symbolic computations, done with the 
help of Maple and our own package Jets [2], are omitted. To stay within the limits given 
by available computing resources we had to restrict the jet order (order of derivatives). 

The answer is given by a third order nonlinear ordinary differential equation (10) 
to govern the functional dependence of the principal curvatures. Incorporation of the 
actual spectral parameter is achieved in Section 3. This can be considered a proof of 
integrability, opening up the possibility to obtain explicit solutions by the methods of 
soliton theory [8, 14, 39]. However, we had to resign ourselves to following this road. 
Neither were we able to establish a Backlund or Darboux transformation [26, 33], which 
would allow us to construct families of exact solutions depending on an arbitrary number 
of parameters. We only remark that seed solutions could be conveniently found among 
the rotational surfaces, see [22, eq. (1)]. 
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The governing equation (10) is explored in Section 4. We identify two basic 
symmetries, scaling and translation (offsetting), and solve equation (10) in terms of 
elliptic integrals. The generic class of integrable Weingarten surfaces we obtained 
depends on one essential parameter (apart from the scaling and offsetting parameters) 
and is believed to be new. In Sect. 5 we establish the integrable Gauss equation (39) 
in the generic case as well as in a number of special cases when the elliptic integrals 
degenerate to elementary functions. All of these special cases could be located in the 
nineteenth century literature. 

Geometrically, surfaces are related by an offsetting symmetry if they are parallel to 
each other, i.e., if they share the same normal line congruence. Therefore, the offsetting 
symmetry indicates that the concept of integrability naturally extends from surfaces to 
their normal line congruences. Section 7 grew out of our attempt to characterize the 
normal congruences of the integrable Weingarten surfaces. We obtain certain relations 
satisfied by suitably chosen metric invariants of the pair of the focal surfaces. Naturally, 
we expect the corresponding focal surfaces to be integrable as well, but a detailed 
investigation had to be postponed to the next paper. 



2. Preliminaries 



We consider surfaces r(x,y), parameterized by the lines of curvature. This is a regular 
parameterization except at umbilic points. The umbilic points are isolated by the 
Hartman-Wintner theorem [18] except for spheres and planes, which are, therefore, 
the only surfaces excluded from consideration. 
The fundamental forms can be written as 
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u 2 dx 2 + v 2 dy 2 , 
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= — dx H dy 2 , 
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(1) 



where p, a are the principal radii of curvature. The radii transform in a very simple way 
under the offsetting symmetry (21) of the integrability problem (unlike the principal 
curvatures p = 1/p, q = 1/a we used in Part I). 

Choosing the orthonormal frame \l/ = (r x /u, r y /v, n), we consider the Gauss- 
Weingarten equations 
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Consequently, the Gauss-Mainardi-Codazzi equations, which are the compatibility 
conditions for (3), read 



v u u 2 v 2 



UUyy + VV XX - ^U X V X - -UyVy + = 0, (4) 



and 



^+ , ap y =o, ^+ f a * =o. (5) 

u p(p — a) v a(a — p) 

As with Part I, we concentrate on Weingarten surfaces, which are characterized 
by the existence of a functional dependence between p and a. We often resort to a 
parametric representation p(w),cr(w) of the dependence. 

Recall that parameters x, y label the lines of curvature; otherwise they are arbitrary. 
In line with Finkel's approach [16], we use this reparameterization freedom to solve the 
Mainardi-Codazzi subsystem (5). The following proposition is a mixture of classical and 
new results. 

Proposition 1. Away from umbilic points, a Weingarten surface can be parameterized 
by the lines of curvature in such a way that 

u = exp / - — - — — dw, v = exp / - — - — — dw. (6) 
J {(T-p)p J {p-o-)a 

The Mainardi-Codazzi subsystem (5) is then identically satisfied, while the remaining 
Gauss equation can be written in the compact form 

Ryy + S xx + T = 0, (7) 

where R,S,T are appropriate functions of the unknown w. Moreover, the constraint 

f---)uv = l (8) 
\p a/ 

can be imposed as an additional condition, and then T = l/(er — p). 

Proof. Writing p(w),a(w) for some function w(x,y), the general solution of the 
Mainardi-Codazzi subsystem (5) is 



u 



= u (x) exp / - — - — —dw, v — v (y) exp / - — - — —dw. 

J {cr-p)p J {p-o-)a 
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Obviously from formulas (1), the multipliers Uq(x), Vo(y) can be removed by an 
appropriate relabelling x — x(x), y — y(y) of the surface's curvature lines. With 
uq = vq = 1, we have 



/ / pa pa \ per 

■uf = exp / I t r — h t — I aw = c , 

./ V p-crp cr -per/ cr-p 



(p - a)p (cr - p) 

where c is an arbitrary constant multiplier. Setting c = 1 by the same relabelling 
argument proves the last relation. 

Having resolved the Mainardi-Codazzi subsystem, we are left with the Gauss 
equation (4) alone. Multiplied by 1/p — 1/c, equation (4) can be written in the compact 
form (7), where 



R 



P 2 , 

— 2 1i dw, 



S = - 



P Jo 
Substituting l/(l/p — l/c) for uv finishes the proof. 



a 2 > ^2 2° ~ P 

—^V dW, 1 — U V 2 2 ■ 



p (J 



(9) 

□ 



3. The classification result 

Employing the Maple package Jets [2], we completed the computer-aided cohomological 
classification outlined in Part I. We have no computer-independent proof of the following 
result. 

Proposition 2. The third-order ordinary differential equation 
in 3 
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determines a unique maximal class of Gauss-Mainardi-Codazzi equations of Weingarten 
surfaces whose initial sl(2, C) -valued zero curvature representation 
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admits a second order formal spectral parameter under the condition that the normal 
form of the zero curvature representation can depend on derivatives of u,v,a,p of no 
higher than the first order. 

Here and in what follows we assume that p is a function of a and the prime refers to 
derivatives with respect to a. A kth order formal parameter A means a power series in 
terms of A up to order k. Part I should be consulted for the other unexplained notions. 

Remark 1. (1) The last proposition provides a complete classification of integrable 
Weingarten surfaces under the following assumptions: The one-parametric zero 
curvature representation takes values in the Lie algebra $1(2), includes the initial zero 
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curvature representation (11) as a member, depends analytically on the parameter, 
and its normal form involves derivatives of no higher than the first order. All these 
limitations can be overcome, in principle [24] , at the cost of requiring significantly more 
computational resources. 

(2) We would like to stress that the only part relying on machine computations is 
the completeness of the classification. All the other proofs in this paper are traditional. 

In the rest of this section we establish integrability of the class determined by 
equation (10). The equation itself will be solved in the next section. 

Proposition 3. The nonremovable spectral parameter exists for all dependences p(a) 
allowed by the governing equation (10). 

Proof. Inspired by the results of the computer-aided classification, we depart from the 
following ansatz for the parameter-dependent zero curvature representation: 



with am, bm, ano, &no, Oi2, «2i, &12 being the unknown functions of a. The problem is to 
solve the zero curvature condition D y A — D X B + [A, B] = for matrix functions A, B of 
u, v, a, p and their derivatives. However, the derivatives are not independent quantities, 
being subject to the Gauss-Mainardi-Codazzi equations. The proper way to deal with 
this situation is to introduce the manifold determined by the equation and its derivatives 
(a difliety [9]). This is fairly easy if the order of derivatives is restricted as it is. Initially 
the derivatives are considered to be independent (jet space coordinates). Considering 
p as a function of a and resolving the Mainardi-Codazzi equations (5) with respect to 
u y ,v x , we can express u y ,v x as functions of u,v,a,a x ,a y . Similarly, the derivatives of 
the Mainardi-Codazzi equations (5) can be resolved with respect to u X y,u yy ,v xx ,v X y, 
giving u X y,u yy ,v xx ,v X y as functions of u,u x ,v,v y ,a,a x ,o- y . Consequently, the Gauss 
equation (4) can be written in terms of u, u x , v, v y , a, a x , a y , a xx , a yy , and then resolved 
with respect to o~ yy . The explicit formulas are somewhat cumbersome, hence omitted. 

With A, B chosen as above, the left-hand side S := D y A — D X B + [A, B] of the zero 
curvature condition S = is a matrix function of u,u x ,v,v y ,o~,o~ x ,o~ y ,o~ xx ,o~ xy . From 
dS/ da xx = and dS/ da xy = we obtain 

bm — frno — a no- 

From either d 2 S/da1 = or d 2 S/do-y = we get a[ u = 0. Hence, a m is a constant, 
which we rename A in anticipation of its role as the spectral parameter. 
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Ap ai2 + a 2 \ , p r ( . . 

OllO = r — r , 6 12 = -7 r 612 + A(a 2 i - ai 2 )] , (12) 

2a{a - p) b 12 a[a - p) 

while <9Sy 9(7^ = can be rewritten as 

a! 12 = 2a 110 a 12 + —t°^ — r(a i2 + 2Xb 12 ), 
Pip - v) 

(13) 

a'21 = — 2anoa2i H — -, \\fl2\ — 2A612), 

P\P ~ <?) 

Modulo these relations, vanishing of S is equivalent to 

612 = —r^ (14) 

We claim that the governing equation (10) arises as the condition that the 
system (12), (13), and (14) be compatible for arbitrary A ^ 0. To prove this, we denote 
P = d\ 2 + CZ21, Q = O12 — 021- 

With ano and b\ 2 taken from formulas (12) and (14), 
respectively, equations (13) turn into 

Pip ~ o) Pip -o) P ip - o) Q 

and the second equation in (12) into 

p\Q 2 - P 2 )Q 2 + p 2 (p' - l)P 2 + 4AV = 0. (16) 

Now the question is whether equations (15) and (16) are compatible. Modulo eq. (15), 
the derivative of (16) with respect to a is 

2p\P 2 - Q 2 )P 2 Q 2 + 2(1 - 3p')p i P 2 Q 2 - 4p'X 2 P 2 

+ (4A 2 + p 2 Q 2 ) [4p'p 2 Q 2 + (p- <j)p" + 2p' 2 - 2p'] = 0. ^ 

This is equivalent to 

[(p _ a )p" - 2p> 2 + 2(1 + 8\ 2 )p']p 2 Q 2 + 4A 2 [(p - a)p" - 2p' 2 - 2p'\ = (18) 

modulo (16), since (18) is the remainder after division of (17) by (16) as polynomials 
in P. Similarly, dividing (16) by (18) as polynomials in Q, we get 

[(p - a)p" - 2p' 2 - 2p'\ [(p - a)p" - 2p' 2 + 2(1 + 8A 2 )p']p 2 P 2 
- 4(1 + 4A 2 )[(a - p) 2 p" 2 - 4p /4 + 8(1 + 8A 2 )p' 3 - 4p /2 ] = 0. 

Differentiating (17) once more and taking the result modulo (15), (19) and (18), we get 
the governing equation (10) immediately. 



Integrable Weingarten surfaces 



8 



Summing up, we obtain a zero curvature representation 
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where P and Q are the square roots to be determined from equations (19) and (18), 
respectively. Away from umbilic points (where p — a), matrices A, B actually exist 
unless (p — a)p" — 2p' 2 — 2p' = when P is undefined. This excludes exactly spheres 
and the linear Weingarten surfaces. The latter surfaces are, however, well known to be 
integrable, being parallel to surfaces of constant curvature (either Gaussian or mean), 
see [38] or [33, §1.5.2]. 

If A = i/2, then we have P = and Q = 1/r 2 , which reproduces the parameterless 
zero curvature representation (11) we started with. □ 

Non-removability of the parameter is ensured by the method [25] (follows from 
nontriviality of the first gauge cohomology group). 

4. Solution of the governing equation 

Apart from the discrete symmetry p -h- a, the governing equation (10) has two obvious 
continuous symmetries, which should be expected in every integrable class of surfaces: 
the scaling symmetry 



p t— > e T p, a t— > e T a 



and the translational symmetry 

pi->p + T, (jhxt + T. 



(20) 



(21) 



The geometric meaning of the latter symmetry is offsetting, also known as taking the 
parallel surface. In terms of position vectors, r is transformed to r + Tn, where n is the 
unit normal vector and T is the distance. 

With the help of these symmetries we can reduce the order of equation (10) by 
two. This can be done by rewriting the equation in terms of the symmetry invariants. 
Since rescaling applies also to the offset, the translational reduction should precede the 
scaling reduction. For the two lowest-order translational invariants we choose 



p-a, 



7] 



(22) 



(recall that the prime denotes the derivative with respect to a). 
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1. If £' = (equivalently, p' — 1), then p — cr = const, which are the surfaces of constant 
astigmatism we dealt with in Part I. 

2. Otherwise, more translational invariants can be computed as derivatives of i] with 
respect to £: 

/ // /// 112 

VP P P 

" S = r^i' ^ = (T^i? " (T^TF ' (23) 

etc. In terms of these invariants, the governing equation (10) reduces to the second-order 
equation 

2^ 2 (v ~ " i\v ~ 3)77| + 2^(r7 - 1)^ " <V + lV = 0. (24) 

As expected, this equation is scaling invariant. To reduce it with respect to scaling, we 
proceed as follows. Besides 77, one more scaling invariant is 

C = C(V - (25) 
Although dispensable, the factor rj — 1 simplifies the computations to follow. 

2.1. If rj' = 0, i.e., p" = 0, then (10) reduces to p' = c, where c is either of —1,0, 1. 
The corresponding surfaces are, respectively, the constant mean curvature surfaces 
(a subclass of linear Weingarten surfaces), the tubular surfaces (surfaces swept by 
spheres of constant radius moving along a space curve) and once more the constant 
astigmatism surfaces. 

2.2. Otherwise p" 7^ and we have 

in 

Cv = y(p~ v) + P' 

In terms of 77, the reduced governing equation (24) becomes the Bernoulli equation 

with the general solution ( 2 = A(rf + 2cq1] + l)i] 2 , where Cq is the integration constant. 
Substituting from eq. (25) yields the separable first-order equation 

^ = ±2 ^L_ vV + 2c oV + 1 (26) 

containing the parameter cq. Being written in terms of the scaling and translation 
invariants, this equation determines the integrable Weingarten surfaces up to rescaling 
and offsetting. Depending on the value of the parameter Co and on the choice of the '±' 
sign, we obtain the following cases. 
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2.2.1 Let c = 1. Equation (26) becomes 

2.2.1.1. With the choice of the plus sign in (27), the general solution is (r? + l) 2 = Ci^ 2 . 
Substituting from eq. (22), we obtain 

(p' + i) 2 = Cl (p-<7)y. 

If ci = 0, the general solution is p + a = const. Otherwise, we apply the transformation 
K = p + a, t = p-o (28) 

to get 

fa«»-4)(!)W 

The equation is separable with a general solution (k — C2) 2 — £ 2 + 4/ci = 0, i.e., 

, 4 I C1C0 
4pa - 2c 2 (p + a) + ^ = 0. 

Ci 

In both cases, c\ = and ci 7^ 0, solutions correspond to the linear Weingarten surfaces. 

2.2.1.2. With the choice of the minus sign in (27), the general solution is (77 + 1) 2 £ 2 = 
C\T). Substituting from eq. (22), we obtain (p' + l) 2 (p — a) 2 = c\p' . For Ci = we have 
the special linear Weingarten surfaces p + a = const again. Otherwise, we apply the 
transformation (28) to get 

d ^ 2 ' 0. 



The solutions are 



k = ± § \^c~i ln(2 v /Z c^C + Vcj - 4c^ 2 ) + c 2 , 
where C2 is the integration constant. 

2.2.1.2.1. For c x < we can write 
V-Ci . 



{ = 2 



or 

p-a 



sinh^± ^ 2 ^ (k - c 2 ) - ln(-ci)^ 
sin h (^ +Co ). 
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2.2.1.2.2. Similarly, solutions corresponding to positive c\ are 

2.2.2 Let c = —1. Equation (26) becomes 

Solutions corresponding to ?y = 1 belong to Case 1 (constant astigmatism surfaces). 

2.2.2.1. The general solution of £,(dr] / d£) = 2i] is rj = c^ 2 . Substituting from eq. (22), 
we obtain the Riccati equation p' = c\{p — a) 2 . 

2.2.2.1.1. For c x > Owe get 

tanh(\/cl a + c 2 ) coth(v / Ci cr + c 2 ) /olX 

p = a t= or p = a -/= (31) 

VCi VCi 



according to whether the integration constant is positive or negative. 

2.2.2.1.2. Similarly, for q < we get 

tan(y/^c^a + c 2 ) co^-/^ cr + c 2 ) , ori . 

p = a or p = a H ^— . (62) 

v-ci v-ci 

2.2.2.2. When solving £(dr} / d£) = —2r/, we get (31) and (32) with p, a interchanged. 

2.2.3. We are left with the generic case c ^ {—1, 1}- Equation (26) has the general 
solution 



G? + c + vV + 2c 7? + l)(c ?7 + 1 + V 7 ^ 2 + 2c r] + l) = c^ ±2 ^. (33) 

If Ci = 0, then 77 = in view of Cq ^ {—1,1}, which yields the tubular surfaces p = const. 
Let us, therefore, assume that c\ 7^ 0. Upon substituting from (22), equation (33) 
becomes a first-order ODE, separable in terms of variables (28) and having the elliptic 
integral 



Vcit ±4 - 2(c + l)(co + 3) Cl t ±2 + (c 2 - l) 2 



as the general solution. The two cases the '±' symbol refers to can be converted one 
into another by the substitution c\ — > (cq — l) 2 /ci. Therefore, we can safely choose the 
sign to be '+', which we do in the sequel. Moreover, if k is a solution, then so is — k (as 
a combination of the p -H- cr switch and a scaling by factor of —1). This is why we often 
ignore the sign of k in what follows. 
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Substituting t — > s/m, m = VJc~i/ (1 — c|)j, we simplify the integral above to 



i _ , , _ , . > /* i±s 2 



m 



/c =-/ ± ( m ^,c), I ± (Z,c)= r = 9 3 ds, (34) 



V 7 ! + 2cs 2 + s 4 



where '±' refers to the signum of Ci/(1 — Cq); in particular, is unrelated to the '±' sign 
in (33). The real parameter c is related to cq by c = ±(c + 3)/(c — 1). 

Formula (34) describes possible dependences p(o~) via the substitution k — p + a, 
£ — p — o. Three independent parameters are involved: m, c and the integration constant 
(the lower limit of the integral). Obviously, m plays the role of the scaling parameter. The 
integration constant can be easily identified with the offsetting parameter T from (21). 

Each dependence between n and £ has a unique representative modulo scaling and 
offsetting, obtainable by fixing the lower limit of the integral 7±(£,c) in (34). This is 
straightforward when c > — 1; we simply redefine i±(£, c) to be 

lis 2 

I±(Z,c)= / I 4 d S . (35) 

Jo Vl + 2cs 2 + s 4 

If, however, c < —1, then the integrand in (34) is real in three separate intervals 
(-oo, -\^u), (- v 7 ^, Vl^), and (v 7 ^ , oo), where 

7± 



= -c ± Vc 2 - 1 > 0. (36) 



We choose the representatives — 7±(— £, c), 7±(£,c), and 7±(£,c), respectively, where 
i±(£, c) is given by (35) in the interval — 7_ < £ < 7-, while 

& c) = / j l± / 2 4 ds, 7+ < £. (37) 
5. Summary of the solutions 

As demonstrated in the preceding section, each integrable class is determined by certain 
relation between the radii of curvature, which can be subject to rescaling p — > Cip, 
er — > ci<7, offsetting p — > p + Co, er — > cr + Co and the twist p -H- cr. 

With the help of Proposition 1, we can find the corresponding integrable Gauss 
equation. To start with, we investigate the generic class determined by formula (34); we 
fix the scaling for simplicity. 

Proposition 4. Assuming 

[* lis 2 

p + a = I±(p-a,c), J±(£,c)=/ , ds, (38) 

J Vl + 2cs 2 + s 4 

the Gauss equation (4) for £ = p — a reads 



R'Zyy + R"Zy + S'U + S"C + T = 0, 



(39) 
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where 

l + ce + A(£,c) c T l 

£ 2 A(£,c) 2 (l + c£ 2 + A(£,c))A(£,c)' 

A(e, c ) = v / i + 2ce 2 +e 4 , T = -\ 

The metric coefficients u,v in (1) are 

2£ ^ ; ' 2£ V 2c ±2 

Proof. We parameterize p and a by £, i.e., we resolve (38) as 

4(e,c) + e /±(e,c)-e 

^ 2 2 

The general form of the Gauss equation, along with the last term T — l/(a — p) = — l/£, 
follow from Proposition 1. To find R',S', we compute 

2 ce 2 +e 4 + ^i + 2ce+e 

i l + 2c£ 2 + £ 4 

2 c^ 2 + e 4 - >/i + 2 C e 2 + e 4 

f l + 2ce 2 + e 4 

from (9) under the constraint (8). These equations need to be integrated once, which is 
easy; the integration constants have been chosen to match equations (8) and (9). Finally, 
from (9) one easily computes the coefficients u, v as u = \ / R'p 2 /p',v= V-S'o-ya'. 

□ 

Apart from the generic class we also obtained a number of special solutions, listed 
in Table 1 (omitting the tubular surfaces). Rows 5b and 6b differ only by translation 
(offsetting) and can be identified one with another. 

The first column contains a determining relation (up to a scaling), while the second 
harbours the corresponding integrable equation in the compact form (7). Table 2 gives 
the principal radii of curvature p, a, metric coefficients u, v, and the variable z (see 
Table 1) in terms of a suitably chosen parameterizing variable w. 

Neither of the special cases is new to differential geometry. Row 1 reflects that, in 
terms of the curvature line coordinates, minimal surfaces correspond to solutions of the 
Liouville equation [5, §351]. Similarly, row 2a reproduces the relation between surfaces of 
negative constant Gaussian curvature and solutions of the elliptic sinh-Gordon equation. 
Row 2b does the same for the hyperbolic sine-Gordon equation and surfaces of positive 
constant Gaussian curvature (or constant mean curvature, by the theorem of Bonnet 
on parallel surfaces). Nowadays, surfaces of constant mean or Gaussian curvature are 
undoubtedly the best understood classes of surfaces integrable in the sense of soliton 
theory (see, e.g., [6, 7, 13, 19, 27, 29] and references therein). 
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relation 








-1 

1. 


p + cr = 









za. 


per = 1 




+ — sinh z = 




OK 
ZD. 


per = —1 


%xx %yy H - sin z 




Od. 


p - o = 


sinh(p + o) 


^Xdllll Z — % )xx \ ^ COX 11 Z — Z jyy -\- CbCll LZ — 


n 
u 


^b 


p-o = 


sin(p + o) 


1 +" q n i —1— 1 pnf y — 1— y i —1— pcp v y — 




4. 


p-o = 


1 


Z xx + i^/ z )yy + 2 = 




5a. 


p-o = 


tanh p 


| (sinh z — z) xx + (coth \ z) yy + coth \ z — 





5b. 


p-o = 


tan p 


| (sin z — z) xx + (cot | z)^ + cot \ z = 




6a. 


p-o = 


cothp 


| (sinh z + z) xa; — (tanh \ z) yy + tanh \ z — 





6b. 


p-o = 


— cot p 


\ (sin z + z) xx + (tan | z) yy + tan | z = 





Table 1. Special integrable cases and the associated integrable Gauss equations 





P 


o 


u 




z 


1. 


w 


—w 


Vw/2 


Vw/2 


— In w 


2a. 
2b. 
3a. 


w 
w 

w + sinh w 


1 

w 

1 

w 

w — sinh w 


w 

vV - 1 

w 

Vw 2 + 1 
w + sinh -U7 


-1 

vV - 1 
l 

\/w 2 + 1 
u> — sinh w 


2 arctanh u> 
2 arctan w 
\w 


2 


2 


2v / coshw — 1 


2 \/ cosh w + 1 


3b. 


w + sin w 


w — sin w 


u> + sin w 


w — sin w 


\w 


2 


2 


2\/l — cos w 


2Vl + cosw 


4. 


w 


w — 1 


w 


(1 - u>)e™ 


e 2w 


5a. 


w 


■u; — tanh w 


w 

sinh w 

w 

sin w 

w 

cosh w 
w 

COS w 


sinh w — w cosh u> 


2w 


5b. 


w 


u> — tan w 


sin w — w cos w 


2w 


6a. 


w 


w — coth w 


cosh W — if sinh w 


2w 


6b. 


w 


w + cot w 


cos u> + w sin if 


2w 



Table 2. Special integrable cases. The radii of curvature p, a, the metric coefficients u,v, and 
the unknown z of the integrable Gauss equation in terms of a variable w. 

It may come as a surprise that the other cases are classical as well. Introduced 
by Weingarten [36, §4] ('eine neue Flachenklasse'), surfaces satisfying the relation 
p — o = sin(p + o) (row 4b) are covered in Darboux [12, §§745, 746, 766, 769, 770] 
('une classe nouvelle de surfaces decouverte par M. Weingarten') and Bianchi [4, §135], 
[5, §245]. Darboux [12, §746] gave a general solution of an equation equivalent to our 
(tanz — z) xx + (cot z + z) yy + esc 2,2 = 0. He also provided a remarkable geometric 
construction in [12, §770], further developed by Bianchi [5, §245]. In a nutshell: the 
middle evolutes are translation surfaces generated by curves of opposite constant nonzero 
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relation limit 



1. 


K = 





J±(£,oo) 


2a. 


K 2 -- 


= £ 2 + 4 


lim 4(m(, 2m 2 )/m 

m=oo 


2b 


K 2 ~ 

IX — 


- -4 


lirn A i (thP — 9?n^ /m 

11111 ± -J- l / / 1*1^ , Cllil fl lit 

m=oo 


3a. 


K = 


arcsinh £ 


lim 1± {mt; , 1 / 2m J /m 

m=0 


3b. 


K = 


arcsin £ 


lim I±(m{;, — \/2m 2 )/m 

m=0 


4. 


e = 


1 


lim I±(m^, — m 2 /2)/m 

m=oo 


5a. 


K = 


— £ + 2 arctanh £ 


7 + (£,-l), |£|<1 


5b. 


K = 


— £ + 2 arctan £ 




6a. 


K = 


— £ + 2 arccoth £ 


iei>i 


6b. 


K = 


— £ — 2 arccot £ 





Table 3. Special integrable cases as limits of i±(£,c) 



torsion; conversely the Weingarten surfaces are orthogonal to the osculation planes 
of the generating curves. Bianchi's research extends to the complementary relation 
p — a = sinh(p + a) (row 3a) as well [5, §246]. The remaining rows (from 4 to 
6b) correspond to involutes of surfaces of constant Gaussian curvature studied by 
Beltrami [3, Ch. 9, §20]. Row 4 (surfaces of constant astigmatism) has been addressed 
in Part I; we have nothing to add except the Beltrami's work as the earliest reference 
we know of. 

Table 3 demonstrates how the cases expressible in terms of elementary functions 
arise as limits of the generic integral (34) for c approaching ±1 or ±oo along a suitable 
curve in the (c, m) space. The tubular surfaces a = const, which are omitted, correspond 
to k — /+(£, 1) = £ + const. 

6. Curvature diagrams 

To exemplify the wealth of classes of integrable surfaces, we plot the representative 
solutions of the governing equation (10) in Figures 1 and 2. We call them curvature 
diagrams, even though the radii of curvature p, a, rather than the curvatures 1/p, 1/er, 
are plotted, contrary to the customary practice [21, Ch. 5]. The benefit is that diagrams 
can be not only scaled arbitrarily, but also freely translated along the dashed line p = a; 
the translation corresponds to offsetting. For clarity, we adjusted the offsetting so that 
the diagrams are symmetric about the origin, i.e., p(a) = —p{—a). 

The diagrams contain plots of functions Za(£, k), Zb(£, k), Ic±(£,k), and 
kl A (£/k,k). All special cases are explicitly included as limits, except the surfaces of 
constant positive curvature (row 2a). These could be obtained as the limit of kI B (£,/k, k) 
as k approaches zero. 

The plots have been calculated using the Legendre normal form [15, 30] of the 
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k = -e~ 2 k = -e~ 1 k = -1 k = -1 




Figure 1. Curvature diagrams n = 2b(£, fc) (the left-hand legend) and n = I\(£,k), |£| < 1 
(the right-hand legend), where k = p + cr, £, = p — o. More can be obtained by rescaling 
and translating along the dashed line p = er, the axis k. Here Xa(£> — 1) = — £ + 2arctan£ 
(row 5b), Xa(£jO) = arcsin£ (row 3b), 2a(£; 1) = 6; Xb(£, — 1) = — £ + 2arctanh£ (row 5a), 
^b(£, 0) = arcsinh£ (row 3a), 2b(£, 1) = £■ Graphs of k = X\(£, &) end on the solid lines |£| = 1. 




Figure 2. Curvature diagrams (a) n — klj^^/k, k), |£| > l/|fc|; (b) /t = 2<7 + (£, fc) (the top left- 
hand legend) and k = —Ic- (£> (the bottom right-hand legend) , where k = p + cr, £ = p — a. 
More can be obtained by rescaling and translating along the dashed line p — a, the axis k. 
In (a), the line k = 1 corresponds to tubular surfaces, k = to surfaces of negative constant 
curvature (row 2b), and k = —1 to the constant astigmatism surfaces (row 4). In (b), 
X c +(£,0) = -£ + 2arctan£ (row 5b) X C -(£, 1) = £ - 2arctan((£ - l)/(£ + 1)) (row 5a after 
reparameterization) . 
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elliptic integrals (35) and (37), which could be of independent interest. As is well 
known, the Legendre normal form depends on the configuration of roots of the quartic 
polynomial II = s 4 + 2cs 2 + 1. 

A) If c < —1, then IT = (s 2 — 7+)(s 2 — 7-) has four real roots Vt± and — Vt± 
given by formula (36). By using the substitution s = y/kr, where k = 7_, we easily 
obtain the Legendre normal form 

, dr, 0<*<1. 

On the right-hand side, we can remove the ± sign from the numerator by allowing k to 
range between —1 and 1. For —1 <£< 1, —1 < k < 1, we have a unified representative 
given by k = Za(£, fc), where 

in terms of the Legendre elliptic integrals E, F. 

For real £ such that |£| > 1, the function Za(£, fc) is complex valued. Yet we obtain 
a real function for l/\k\ < £ by choosing the lower limit of the integral to be 1/k, 
-1< fc < 1. Thus, 



ii/w V(i-r 2 )(i-fcV) v fI ' |fc| 

-z A (-e,fc), e< ~wi 



B) Similarly, when c > 1, then 7 ± < 0, the roots V7±, — V7± of II are purely 
imaginary, and 

1 r /V /T ^ + ^ ii^ 2 , n i i 

Vt ( ~2*~) = i, TcTTTjirTlvT 0<fc<h 

The two representatives can be unified into « = Zb(£, fc), where 

Zb&*)=/ 77- = ; | 2 dr = -g(gi; fc) + — F(£i; fc) 
Jo V(l +r 2 )(l + fc 2 r 2 ) fci k\ 

for -1 < fc < 1. 

C) When — 1 < c < 1 (four distinct complex roots), we substituted 

l+Vkr 

s = T37fc? 0<k<1 > 
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to obtain two more representatives k = Xc+ an d k = Zc_, where 



fJc±(t,k)- J c ±(0,fc), £>0, 
\-2b±(-e,fc), e<0, 



V ^ ; 1 + e (A; + l)i \£ + lVk ' y i V^ + lv^ 

(1 =b l)fc - 3 =b 1 fjfe-1, 

£± = 



2 I 2. 

A; 2 - 6/c + 1 
(^ + 1) 2 ' 



7. Normal congruences and their focal surfaces 

The fact that the governing equation (10) has the offsetting symmetry (21) is not a 
pure coincidence. Being invertible, the offsetting transformation r t— > r + Tn preserves 
integrability in every reasonable sense of the word. Surfaces related by the offsetting 
transformation are said to be parallel and either all are integrable or none is. However, 
parallel surfaces can be alternatively described as normal surfaces to the same line 
congruence. Consequently, integrability is a property of this congruence and, therefore, 
must have an expression in terms of congruence invariants. 

Normal congruences of Weingarten surfaces, also known as W^-congruences, are 
rather special with regard to properties of their focal surfaces. It is therefore natural to 
look for characterization of the former in terms of the latter. Naturally, we expect the 
focal surfaces of integrable W^-congruences to be integrable as well. 

Recall that a generic surface has two focal surfaces (often considered as two sheets 
of a single surface), 

i-W = r + en, = r + pn. 

each of which is formed by the evolutes of one family of the curvature lines. Focal surfaces 
can degenerate into a line or even a point. In the case of a Weingarten surface r with 
fundamental forms (1), one of the focal surfaces degenerates into a line if a y = a'w y = 
or p x = p'w x = 0; both degenerate into a point if the surface is a sphere (already 
excluded from consideration); otherwise they are regular surfaces. Therefore, we assume 
p'a' in what follows. 

To compute the respective first and second fundamental forms 1^ and IlW, % — 1, 2, 
we proceed as follows. In view of the Gauss-Mainardi-Codazzi equations (4) and (5), 
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the Gauss-Weingarten (3) equations can be written as 
u x ap'u 2 w y u 2 1 

U p{p — <7)V p p 

ap'w y pa'w x 

xy p(a-p) x + a(p-a) y ' 1 } 

pav 2 w x v y v 2 1 

yy a (a - p)u 2 v y a y a y 

One easily finds 



= - — -r x + a'w x ti, rft = a'w v ti, = — 

p y v 

r {2) = p'w x n, r y 2) = - — -r y + p'w y n, n (2) = — . 

(J u 



Using the equations (40) and (1), we get 

IW = ip -f u2 dx 2 + d. 2 , = dp 2 + ip -f v2 dy 2 , 

p a 



(41) 



where dp = p' dw = p'(w x dx + w y dy), da = a' dw = <j'(w x dx + w y dy). 
With u,v determined from Proposition 1, we can write 

I« = (fU(a) dx) 2 + da 2 , l& = U {2 \P) dy) 2 + dp 2 . 

Hence, all focal surfaces rW corresponding to a given dependence p(a) are isometric. 
Moreover, the first fundamental forms (41) are typical of surfaces of revolution. These 
are among the classical results by Weingarten [36]. 

Omitting details, we further compute the second fundamental forms 



II^^f^d^-^dA ll^ = - P -^(^dx 2 -^dy 2 ) (42) 
v \ p a / u \ p a J 



and note that they are conformally related, which is another way to express Ribaucour's 
classical result [32] that asymptotic coordinates on and r^ 2 ^ correspond. The 
Gaussian curvatures are 

" detlW " (P-^)V' A ~ detl< a > ~ (p-a)V [46) 

Consequently, the focal surfaces have one and the same sign of the Gaussian curvature, 
which we denote as e. We have e — — 1 (both focal surfaces are hyperbolic) if and only 
if dp/ da = p'/a' > (if p increases as a increases), and +1 if dp /da < 0. The relation 

KmKm - i^hr (44) 
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away of umbilic points is known as the Halphen theorem (see [4, §129]). 

As we have already explained, to every particular relation p(o~) of curvatures there 
corresponds an isometry class of focal surfaces, which contains a unique rotational 
representative (which is the way the classes have been characterized in the classical 
literature). However, we believe that a description in terms of metric invariants is more 
appropriate. It is convenient to choose 

«« = ' 



where eK® = \K^\ is the absolute value of the Gaussian curvature of the ith focal 
surface. 

Further, let 7W be defined by 

{l) _ {p-a){p"a'-a"p')-2p'a\p'-a') 
7 ~ 2{-ep'a'f' 2 

(2 ) (p-a)(p"a'-a"p') + 2p'a'(p'-a') 



2 (-ep'a'Y' 2 



One can directly check that \^\ equals the norm of the gradient of «W with respect 
toI«, 

| 7 (i) | = ||grad (<) K (<) || (i) = \// {l) (grad (l) K «,grad (<) k w ) . 

Hence, 7W is a metric invariant of the respective focal surface. It is sometimes more 
convenient to use invariants 

_ [( P -*)(p"*'-*"p')-2p'aV-a')] 2 
16 {p'a'f 

' (46) 
r{2) _ [{p-a){p"a'-a"p') + 2p'a'{p'-a')f 

16 {p'a'f " " ' 



satisfying 

7 W2 = _4 e G(0 > _ l6G w ^ (i)3 = 7«(grad w grad (i) fsT (i) ). 

Clearly, both «« and G« are functions of w. Consequently, can be considered as 
a function of unless «W is a constant. Our nearest aim is to establish the dependence 
between k^> and in terms of the dependence between p and a. 

Proposition 5. Let the principal radii of curvature p,a of an integrable surface satisfy 
the generic relation (34). Then the metric invariants and satisfy the relations 

G «=ii ± ^:^)(-i + ^:i K « 

V c =F 1 m / \ V c =F 1 m 
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Furthermore, 

gUgw = (—) 2 

is constant (hence, so is the product ). 

Table 5 lists the product G^G^ and the algebraic relations between G^ and 
in the special cases. 

Proof. For simplicity, we start assuming a fixed scaling, i.e., we depart from formula (38). 
We routinely compute 



2 , . A , r, .2 I 4^2 



(1) {l±w z + Vl + 2cw z + wV (2) U±w"- Vl + 2c^ + w% 

2(c T l) W 4 ' 2(ct1)«j 4 

Consequently, e = sgn(c =F 1)), and 

(1) _ 1 ± w 2 - \/\ + lew 2 + w 4 {2) _l±w 2 + \T\ + lew 2 + w 4 
a/2 |c =f 1 1 ' ^ v^2|c=f 1| 

Furthermore, 

{1) _ (lTw 2 + y / l + 2cw 2 + w 4 ) 2 (2 ) (l Tw 2 - y / l + 2cw 2 + w 4 ) 2 
2(c T lV 2 ' 2(cTl)w 2 

Under the scaling by factor of m, the metric invariants and become /m 2 and 
mn' 1 ', respectively, while G^ remains invariant. Formulas (47) are then easily checked. 
Moreover, all three metric invariants are invariant under the offsetting (21). 

Formulas for and k;W in the special cases are given in Table 4 along with the 
sign e of the Gaussian curvatures. □ 

Summarizing, focal surfaces of integrable Weingarten surfaces belong to the 
isometry classes specified in Proposition 5. 

A natural question is whether the condition G^G^ = const or, equivalently, 
syW^i 2 ) — const, is not only necessary, but also sufficient for the condition (10) to 
hold. 

Proposition 6. Under the condition 7^ + 7^ ^ 0, a surface satisfies the governing 
equation (10) if and only if the product 

7 (D 7 (2) = ± || grad (l) K (l)||(l) || grad (2) K (2)||(2) (4g) 

is constant. 

Proof. Assuming the p(a) dependence, 7^ +7^ simplifies to (p — a)p"/ v^p'] 3 and the 
product in question to 

' 7 ep' isp' 3 ■ 
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(1) 



K 



(2) 



G™ 



1. 

2a. 

2b. 

3a. 

3b. 

4. 

5a. 

5b. 

6a. 

6b. 



2 \w\ 
1 

~2 ~~ 



+ 1 



w 

-1 + coshw 

1 — cos w 
1 

tanh 2 w 
tan 2 iu 
coth 2 w 



cot 2 w 



2 H 
\w 2 - 1| 

w 2 + 1 

1 + cosh w 

1 + cos w 
1 
1 

1 

1 

1 



-1 
1 

l 

1 + cosh w 

1 — cosh w 
1 + cosw 

1 — cos w 

1 

sinh w cosh u> 



sin 2 w cos 2 u> 



sinh 2 w cosh 2 w 
1 

sin 2 w cos 2 u> 



-1 

- W 2 

w 2 

1 — cosh w 

1 + cosh iu 
1 — cos w 

1 + cos to 









Table 4. Special integrable cases. Metric invariants of focal surfaces in terms of w. 



e G^G™ gU( K W) G {2 \^) 

T I -i -l -i 

2a. 1 1 -l±/t (1) -1±k {2) 

2b. -1 1 -1 + kW -lW 2) 

2 2 

3b. -1 1 -l + 4n -1 + 



4. -1 



5a. -1 ( r =) 

5b. 1 -fv^ + ^^l 

6a. -1 f V 7 ^ i=] 

6b. 1 

v v^y 

Table 5. Special integrable cases. Relations between metric invariants of focal surfaces. 
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Factorizing the a-derivative of this expression as 

2ep V 2p p — (j \p — a) / 

and comparing to the governing equation (10) proves the proposition. □ 

It follows from the proof that condition (48) also holds when p" = 0, i.e., if there 
is a linear relation between the radii of curvature. As of now, there seems to be no 
indication towards integrability of the latter class (except when p ± a — const, which 
satisfies (10) as well). 



8. Conclusions and future work 



In this work we singled out a class of Weingarten surfaces on the basis of its solitonic 
integrability. Although special cases were not unknown to nineteenth century geometers, 
the overall result appears to be new. We also characterized integrability in terms of 
metric invariants of the focal surfaces. 

For time reasons, many questions had to be left for further research. We do not know 
the Backlund transformation, recursion operator, bi-Hamiltonian structure and other 
attributes of integrability. We did not provide any solutions to the Gauss equation (39). 
We do not know what is the true geometric meaning of the spectral parameter. Even 
the task of computing third order symmetries of the Gauss equation proved to be too 
complex. 

We have seen in Part I that integrability of surfaces of constant astigmatism is 
attributable to the fact that their focal surfaces are pseudospherical. In the general 
case, the existence of an integrability-preserving relation to previously known integrable 
surfaces is an open problem. 

Our nearest goals include exploring the induced Bianchi type transformation 
between surfaces satisfying relations (47) as well as investigating the extended 
symmetries of the class in the sense of Cieslihski [10, 11]. 
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